Field theory of the inverse cascade in two-dimensional turbulence 
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A two-dimensional fluid, stirred at high wavenumbers and damped by both viscosity and linear 
friction, is modeled by a statistical field theory. The fluid's long-distance behavior is studied using 
renormalization-group (RG) methods, as begun by Forster, Nelson, and Stephen [Phys. Rev. A 16, 
732 (1977)]. With friction, which dissipates energy at low wavenumbers, one expects a stationary 
inverse energy cascade for strong enough stirring. While such developed turbulence is beyond the 
quantitative reach of perturbation theory, a combination of exact and perturbative results suggests 
a coherent picture of the inverse cascade. The zero-friction fluctuation-dissipation theorem (FDT) 
is derived from a generalized time-reversal symmetry and implies zero anomalous dimension for the 
velocity even when friction is present. Thus the Kolmogorov scaling of the inverse cascade cannot 
be explained by any RG fixed point. The (3 function for the dimensionless coupling g is computed 
through two loops; the g 3 term is positive, as already known, but the g 5 term is negative. An ideal 
cascade requires a linear (3 function for large g, consistent with a Pade approximant to the Borel 
transform. The conjecture that the Kolmogorov spectrum arises from an RG flow through large g is 
compatible with other results, but the accurate k~ 5 ^ 3 scaling is not explained and the Kolmogorov 
constant is not estimated. The lack of scale invariance should produce intermittency in high-order 
structure functions, as observed in some but not all numerical simulations of the inverse cascade. 
When analogous RG methods are applied to the one-dimensional Burgers equation using an FDT- 
preserving dimensional continuation, equipartition is obtained instead of a cascade — in agreement 
with simulations. 

PACS numbers: 47.27.Gs 



I. INTRODUCTION 



The cascade of energy to low wavenumbers in two- 
dimensional turbulence [lj, more than other turbulence 
problems, is suited to the standard methods of statis- 
tical field theoryand the renormalization group (RG). 
These methods Q, originating in quantum field theory, 
show that arbitrary short-distance interactions lead to 
long-distance behavior described by a local, renormaliz- 
able action — an effective field theory. Correlation func- 
tions computed from this action contain ultraviolet (UV) 
divergences that can be eliminated by redefining the pa- 
rameters and fields. The divergences leave their mark, 
however, in the dependence on the renormalization scale 
and the resulting anomalous scaling laws. 

The classic application of statistical field theory is to 
critical phenomena (second-order phase transitions) in 
condensed matter [2j. The infrared (IR) scale invariance 
of correlation functions at the transition temperature is 
explained by a fixed point of the RG flow. The inverse 
energy cascade of two-dimensional turbulence is likewise 
believed to be nearly scale invariant, and one might sus- 
pect that a similar fixed point is responsible. We will see, 
however, that no fixed point can reproduce the observed 
£.-5/3 ener gy spectrum. Rather, we will argue that the 
inverse cascade arises from a nontrivial RG flow and thus 
is not expected to be completely scale invariant. 

In the study of turbulence, a deviation from scale in- 
variance (a dependence of dimensionless physical quan- 
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tities on scale) is referred to as intermittency Q. While 
intermittency is recognized as a property of the three- 
dimensional direct cascade of energy, its existence in the 
two-dimensional inverse cascade is unsettled. For a non- 
stationary inverse cascade, in which energy is not dissi- 
pated but progresses to ever-lower wavenumbers, inter- 
mittency is not observed in numerical simulations 0, ; 
a theoretical explanation has been given |fj. In this pa- 
per we deal solely with the stationary regime, where an 
inverse cascade requires a low-wavenumber energy sink. 
With few exceptions, simulations of such a cascade con- 
firm the fc -5 / 3 energy spectrum initially predicted [lj on 
the basis of scale invariance. But one set of simulations 
7] finds strong intermittency in fourth- and higher-order 
velocity correlations. Other simulations ||| and experi- 
ments 9], though, find no significant intermittency. The 
various studies differ mainly in the precise form of the 
dissipation terms. The evidence suggests that intermit- 
tency in the stationary inverse cascade, permitted by our 
theory, is at least possibly realized. 

It is our restricted focus on the inverse cascade that 
allows us to work with a purely local field theory. The 
random force that stirs the fluid is correlated over a lim- 
ited range and is effectively local in a long-distance de- 
scription. As with the RG treatment of quantum fields 
and condensed matter, we expect all short-distance de- 
tails to become irrelevant except as they are manifested 
in local, renormalizable couplings. The two-dimensional 
direct cascade of enstrophy to high wavenumbers [| thus 
falls outside our scope. A previous RG analysis of two- 
dimensional turbulence [lOj is formally similar to ours, 
but it follows three-dimensional studies by adopting the 
long-range force correlations necessary for a direct cas- 
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cade; even its derivation of the inverse cascade relies on 
nonlocal forcing. Here we apply RG methods in the fa- 
miliar domain of local field theory, which should allow a 
physical treatment of the inverse cascade. The explana- 
tion of the direct enstrophy cascade may rest on entirely 
different foundations, such as conformal invariance 

The work most closely aligned with our theoretical ap- 
proach is due to Forster, Nelson, and Stephen (FNS) [l2fl . 
At the technical level, our contribution is to add linear 
friction to FNS model A in d = 2 and compute the RG 
flow to the next order of perturbation theory, two loops. 
The inclusion of friction, which dissipates energy at low 
wavenumbers, makes our theory capable in principle of 
describing the inverse cascade and its fc -5 / 3 spectrum — 
unlike FNS model A, which gives a k 1 spectrum corre- 
sponding to energy equipartition in d = 2. We also note 
a difference in our viewpoint from that of FNS and oth- 
ers 0, 0, 0] who apply RG methods to turbulence by 
analogy with critical phenomena. These authors seek a 
controlled IR-stable fixed point by starting with a loga- 
rithmically divergent field theory and then decreasing the 
dimension of space or the exponent of the stirring-force 
correlation by e. This adds to the (3 function a negative 
linear term proportional to e. With the usual positive 
one-loop term, there exists an IR-stable fixed point at a 
coupling that goes to zero with e; the fixed-point theory 
can then be expanded in e instead of the original cou- 
pling. Like FNS, we work in d = 2 — e to regulate UV 
divergences, but we ultimately take e = 0, so that the 
fixed point is trivial. Our inverse-cascade model lies not 
at a fixed point but in the region of large dimensionlcss 
coupling. 

Naturally the use of perturbation theory is question- 
able for strong coupling. Our perturbative results will 
have direct quantitative application only to the extreme 
IR limit controlled by the trivial fixed point, which is of 
some interest in itself. Nevertheless, we will make reason- 
able conjectures about the theory's strong-coupling be- 
havior that are consistent with the inverse cascade, bear- 
ing in mind the dangers of the nonperturbative regime. 
Besides the concern with the numerical accuracy of ex- 
trapolations, there are fundamental difficulties at strong 
coupling. The anomalous dimensions of operators may be 
large, and the relevance of terms in the action may dif- 
fer from the weak-coupling case. Furthermore, at strong 
coupling, there is no simple relation between the cou- 
plings in different renormalization schemes, such as the 
Wilsonian cutoff (useful for physical interpretation) and 
minimal subtraction (convenient for systematic calcula- 
tion). We may hope that these subtleties do not affect 
the main conclusions even at very strong coupling. At 
least we know that the theory of critical phenomena in 
d = 4 — e is extrapolated to e = 1 (moderate coupling) 
with acceptable results. 

In Sec. [H] we describe the basis of our theory and our 
method of calculation, confirming the one-loop RG flow 
of FNS |2 . In Sec. IIIII we present symmetries and other 
properties of the theory that do not involve a dubious 



extrapolation to strong coupling. In Sec. lIVI we compute 
the two-loop term of the (3 function. In Sec. [V] we relate 
the plausible strong-coupling behavior of the theory to 
the phenomenology of the inverse cascade. In Sec. IVII 
as a test of our methods, we consider a rather different 
model, the UV-stirred one-dimensional Burgers equation. 
A summary and discussion are presented in Sec. IVIII 

II. FRAMEWORK 

A. Path integral for the Navier— Stokes equation 

The Navier-Stokes equation for the velocity field Vj of 
an incompressible two-dimensional fluid is 

Vj + v k V k Vj + VjP - vV 2 Vj + aVj = fj, (1) 

where P is the pressure divided by the density, v is the 
kinematic viscosity, a is the friction coefficient, and fj is 
the force per unit mass. The incompressibility condition 
=0 allows the velocity to be expressed as 

Vi = eijVjip, (2) 

where ip is a pseudoscalar field called the stream function 
and Cij is the alternating tensor, which in two dimensions 
satisfies 

fy-ffcl = SikSjl — SuSjk- (3) 

Upon writing Eq. Q in terms of i\) and applying the 
operator e^-V,, we obtain [T(il | 

-V 2 tp - eijViV^Vjtp + vV A i) - aV 2 ip = 17. (4) 
Here 

V = V x f , (5) 

with the notation 

a x b = tijdibj. (6) 

In the real three-dimensional world, Eq. Q is a good 
approximation for a thin fluid film provided either (a) 
the boundary surface(s) and the coordinate system are 
rotating rapidly about a perpendicular axis or (b) the 
fluid is conduc ting and subject to a strong perpendicular 
magnetic field In either case, boundary-layer effects 
produce a linear friction parametrized by a, which has 
units of frequency. 

For convenience, our field-theory calculations will use 
the method of dimensional regularization based on 
continuation to a noninteger spatial dimension d = 2 — 
e. Because in the end we are concerned only with d = 
2, we adopt a formal continuation of the Navier-Stokes 
equation that preserves its two-dimensional features. For 
general d, we retain the stream-function representation 
by choosing a "physical" two-dimensional subspace that 
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contains the tensor and the external wavevectors of 
correlation functions. Denoting by 0^ the projector onto 
the physical subspace, we now have 



(7) 



We take Eq. J3J as the equation of motion for i/>, with 
V 2 = VfcVfc interpreted as the c£-dimensional Laplacian. 
In this way we preserve (for v = a = r\ = 0) the formal 
conservation of the energy and enstrophy, 



Q 



d d xV 2 tpV 2 ip. 



(8) 
(9) 



We assume that the fluid is stirred by a Gaussian ran- 
dom force that is uncorrelated in time |T^ . with 

(/ i (w,k)/ J (w',k')) = (27T) d+1 6(u, + u,')S(k + k') 

x {6 lJ -k l k :j /k 2 )D(k 2 ), (10) 

so that r\ is also Gaussian with 

(rj{u, k) 7](lu', k')> = (2ir) d+1 6(uj + uf) 6(k + k') k 2 D(k 2 ). 

(11) 

A classical system with random forcin g ca n be treated in 
the formalism of quantum field theo ry |17| , including the 
path-integral representation [Hi, ll9T |20|. Upon introduc- 
tion of a pseudoscalar field p conjugate to rj, correlation 
functions of are given by the path integral 



with the action 



Vi/iT>pF[ip] e~ 



(12) 



S= / dtd d x [\{-V 2 p)D(-\7 2 )y 



= J dtd d x[\{~V 2 p)D{-\7 2 )p + iup\/ i ^ 

- iap\J 2 i> - ipW 2 ip - ieijViVkpVjipVkif] , (13) 

where we have integrated the pipip term by parts. The 
Jacobian determinant from changing variables from 77 t o 
if> is an unimportant constant by virtue of causality [2fJ • 

B. Relevance of couplings 

In the field theory based on the action (|13|) . the long- 
distance behavior is governed by just the renormalizable 
terms — those with coefficients whose scaling dimensions 
with respect to wavenumber in d = 2 are zero (marginal) 
or positive (relevant) 0- We assign scaling dimensions 
d^j, d p , d t to the fields ip and p and to the time t by 
requiring that the highest-derivative quadratic terms in 
the action have dimensionless coefficients in d — 2, since 



these terms control the asymptotic behavior of propaga- 
tors and thus the UV convergence or divergence of dia- 
grams. Because of the additional derivatives, the viscous 
term ivpV^ip 1S clearly less relevant than the friction term 
— iap~S/ 2 ip. In fact, it is commonly said that viscosity is 
irrelevant to the inverse cascade, but we now show that 
this cannot be taken in the technical sense. If the viscous 
term is ignored, then the pS7 2 ip and pV 2 V> terms have di- 
mensionless coefficients only if dt = and d^p = —d p . 
For a nontrivial theory, the ptptp term must be renor- 
malizable, giving 2 > 4 + d p + 2d^, = 4 — d p . With 
dp > 2, there exists no local renormalizable forcing term 
quadratic in p. 

Let us therefore retain the viscous term and recom- 
pute the scaling dimensions. The p^7 4 ip and pV 2 ip terms 
possess dimensionless coefficients only if dt = —2 and 
d^p = —dp. Renormalizability of the p^pip term now gives 
4 > 4 + dp + 2d^ = 4 — dp, so that d p > 0. It remains 
to specify the forcing term. We assume that the exter- 
nal forcing is confined to a band of high wavenumbers, 
as in model C of FNS (l2^|. The effective forcing at low 
wavenumbers is generated by renormalization; because 
the interaction in Eq. i|13|) has two spatial derivatives act- 
ing on p, at least two derivatives must accompany each 
factor of p in any term so generated. The only renormal- 
izable forcing term is then \D§\I 2 p\I 2 p, whose coefficient 
is dimensionless for d p = d^p = 0. This effective forcing 
has D(k 2 ) = D k 2 , as in model A of FNS. All terms in 
the action are now marginal, except for the friction term, 
which is relevant (coefficient of dimension 2). The only 
other renormalizable terms that could be generated are 
ones containing only tj), but these are not generated (see 

sec, urn. 

Next we label the fields and the time in Eq. 1)1 3[) with 
the subscript "phys" and introduce rescaled variables to 
simplify the action. Tentatively seeking to set all coeffi- 
cients other than forcing and friction equal to 1, we take 

V'phys = MP, Pphys = V~ X p, t phys = V~H. (14) 

The result is 



S = J dt d d x (y 2 V 2 P V 2 p + i P V 4 ^ 



iv 1 apV 2 il> - ipV 2 %l> - iV k pS7 jipV ki>) , 



where 



(15) 



(16) 



We adopt, however, a different rescaling that will be par- 
ticularly convenient in Sec. IIII Bl 



The final form of the action is then 



(17) 



S = J dtd d x (y- 1 X7 2 p\7 2 p + ig^pX/ 4 ^ 

- ippV 2 ip - i P V 2 Tjj - ieij ViVfcpVjV'VfcV'), (18) 
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k 2 (co + ip+ig- l k 2 ) 
-l 



(WW) 



— ~T 



+ (p+g~ l k 2 ) 2 




-r w = /(^-^) kl xk 2 



FIG. 1: Feynman rules for the action I18|l . 



co, k — q 
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where 



co, k 



p= (g^a = l7 V /2 a 



(19) 



For general d = 2 — e, the scaling dimensions in Eq. p8|) 
are 



— dp — dt — ~ 2 + ie, 



(20) 



In two dimensions, g is a dimensionless coupling and p is 
analogous to a mass parameter in quantum field theory. 



C. Feynman rules 

Correlation functions can be calculated for the action 
(jl8|l using Feynman diagrams whose lines carry both fre- 
quencies and wavevectors. We represent the fields tp and 
p by wiggly and plain lines respectively. The ingredients 
of the diagrams, shown in Fig. ^ are the propagators 
(pip)o, (ipip)o, obtained from the quadratic terms in the 
action, and the vertex factor — TX^ P , obtained from the 
cubic term. We label these quantities with the subscript 
because they are the tree-level contributions to the ex- 
act two-point correlation functions (pip), (ipip) an< i the 
exact three-point one-particle-irreducible (1PI) function 

The remaining contributions arise from diagrams con- 
taining loops. For these diagrams, we integrate over each 
loop frequency and wavevector according to 



duj d d k 
2^ (2n) d 



(21) 



Because the integrand is a rational function of the fre- 
quencies, the u) integrations can easily be performed by 
the contour method before integrating over wavevectors. 
This contour integration shows that a 1PI diagram (or 
sub-diagram) vanishes if all its external lines are ip |20j| - 
since there is a closed loop of (pip) a propagators and the 
integrand is an analytic function of the loop frequency in 
the upper half-plane. 



FIG. 2: One-loop diagrams for the two-point 1PI functions 
— r^ p and — r pp at zero frequency. 



Unlike many field theories in a low number of spatial 
dimensions, ours does not contain IR divergences even 
for p = 0. This is because, after integration over frequen- 
cies, internal-line propagators scale as fc -2 , but at least 
one further factor of the wavevector arises from each of 
the two vertices that a line connects. Hence the inte- 
grand does not diverge as the wavevector of any internal 
line goes to zero. The frequencies and wavevectors of the 
external lines act as an IR cutoff. For simplicity, our cal- 
culations will adopt another IR cutoff by assuming that 
p > 0; then 1PI diagrams are analytic at zero external 
frequencies and wavevectors. 



D. One- loop renormalization 

The coefficients of the quadratic terms in the action 
(with frequency \ an< i wavevector q) are 



ri p = q 2 ( x + ip + l g- 1 q 2 ), 



(22) 
(23) 



These are corrected at one loop by the two-point 1PI 
diagrams in Fig. [21 Because of the external wavevectors 
in the vertex factors, the diagrams are 0(q i ), and so we 
can set x — 0- We now demonstrate the computation of 



the — rf p diagram, to show the basic methods to be used 



for two-loop diagrams in Sec. IIVI 

The frequency integral of the propagators is 



I : 



(111; 
2^ 



[g- l k- 2 (uj + ip + ig^k 2 )- 1 



x(Lj + ip + ig X |k-q| 2 ) 1 (io-ip-ig X |k-q| 2 ) X ] 
~ig/Ak 2 



(gp + k 2 - k • q + \q 2 )(gp + k 2 - 2k- q + q 2 )' 



(24) 
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obtained conveniently by closing the contour in the upper 
half-plane and picking up one pole. We next multiply by 
the vertex factors and expand to 0(q 4 ): 

i(k 2 - 2k • q)(k x q) j(2k • q - q 2 )(k X q) I 

= (k 2 - 2k • q)(q 2 - 2k ■ q)[k 2 q 2 - (k • q) 2 ]/ 
6i<?(k-q) 2 [fc 2 g 2 -(k-q) 2 ] 



A(gp + k 2 ) 3 
iff[fcV+4(k.q) 2 ][fc 2 g 2 -(k-q) 2 ] 
Ak 2 (gp + k 2 ) 2 



, (25) 



where k 2 = Gijkikj is the squared projection of k onto 
the physical subspace and we infer from Eq. J7J that 



(k x q) 2 = k\q 2 - (k • q) 2 



(26) 



We omit 0(q 3 ) terms in Eq. I)25|l because they will now 
disappear when we average over directions of k. 
With the d-dimensional isotropization formulas 



ft/ j, ftj q ' r f\ — ^ ft/^ f\ j tljfc ft<£ ' r f\i 



which imply 



(k.q) 2 Af 



2k 2 

T' 



(k-q) 
r, (k-q) 4 



d(d + 2) 



2 k 2 q 2 



d ' 
3fcV 



(27) 



(28) 



d(d + 2)' v " My ' d(d + 2)' 
the value of the diagram to 0(q A ) becomes 



-rf = -igq 4 



dk 2TT d / 2 k d - 1 fc2 (6 + + dfc 2 



o (2tt)<* r(|d) 4d(d + 2)( 5P + fc2)3 



6-d 



-igq 



(2 + d)2 4+d r(id)sin(i7rci) V#P 



l-d/2 



(29) 



For d — 2 — e with e — > 0, we find 



-if = igq 4 



1 ln( 5 p/47r) + 7e - 1 



327T6 



64tt 



+ 0(e) , 



(30) 

where 7e is the Euler constant. The result for the other 
one-loop diagram (taking into account the symmetry fac- 
tor) is very similar: 



_ t pp 



gq 



l 



32tt£ 



ln(ffp/47r) + 7e 
64tt 



+ 0(e)). (31) 



We have carefully obtained the 0(e°) terms, which will 
be important in Sec. IIVI 

The method of minimal subtraction plj expresses the 
bare couplings g and p. of respective dimensions and 



2 — he, in terms of a renormalization scale p. (dimension 
1) and dimensionless renormalizcd couplings g and p as 



9 

9 : 



■I 2 



[9- 



g 3 a!6 1 



o(f)}., 



(32) 
(33) 



where the corrections (called counterterms) involve only 
negative powers of e. The coefficient of the 0(g 3 ) coun- 
terterm, and the absence of any counterterms for p, are 
obtained by requiring that T^ p and T pp to 0(q 4 ) be finite 
at e = when expressed in terms of g and p: 

T^ = q 2 ( X + ip) 



iq 



327T6 



+ 0(£ u ) +0(g 3 ) 



= q 2 (x + ^^ /2 -p) 



iq 4 p'^ 2 



*k - a i+o^)) + om 



62ne e 



Hence 



Oi 



1 

32^' 



(34) 



(35) 



(36) 



and no counterterms of any order are needed for p be- 
cause the q 2 term of T^ p is not renormalized. 

As for the three-point 1PI function -T^ p ', we will 
see in Sec. IIII A"l that it is related by Galilean invariance 
to the q 2 x term of r^ p . Because the latter term is not 
renormalized, we have 



v 4>4>p 
_i o 



i{kl 



fc|)ki x k 2 , 



(37) 



up to irrelevant terms with more factors of wavevector. 
We have thus rendered the theory finite at one loop by 
renormalizing only the coupling g, without the need for 
counterterms to rescale the fields ip an d V or the time 
t. Crucial for this was the equality of the l/327re terms 
in Eqs. I|3U|I and l|31|l . We conclude that the anomalous 
dimensions are zero at one loop, and in Sec. IIII 51 we will 
show that in minimal subtraction they are exactly zero, 



7V = Ip = It = 0. 



(38) 



In a different renormalization scheme, or with a different 
definition of the fields and the time such as Eq. (|14|) , the 
anomalous dimensions would not vanish identically, but 
at any RG fixed point their values are universal Q and 
so they would be zero there. 

The (3 functions for the dimensionless couplings are 
determined by the RG invariance of the bare couplings, 



= 0. 



(39) 
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We obtain 

P(d)=-h9 + ^- + 0(g 5 ), (40) 
/3(p) = (-2 + £e)/5. (41) 

The force correlation adopted by FNS ^1 differs from 
Eq. IjlOl) by a factor of 2, and consequently the dimen- 
sionless coupling of FNS is A = 2~ 1 / 2 g. Thus we have 
confirmed the FNS result 

/3(A) = -I e A+^- + 0(A 5 ). (42) 
III. GENERAL PROPERTIES 



the forcing, viscous, or friction terms, which contain ip 
only as W 2 ip an d are trivially invariant. We are left with 

SS = -Jdtd« x (ipu^W + fc^Vtf^V^ 

+ ie i:j \7iS/ k pS7 ji) e kl u{). (46) 

Upon integration by parts and use of Quui — itj, the first 
two terms cancel and the third vanishes. 

In Fourier space, the transformation (|45|l becomes 

5tp(uj, k) = u • k ^-ip(uj, k) - i(2ir) d+1 d(Lu) u x V<5(k), 

UU) 

ip(w,k) = u.k| ; p(w ) k). (47) 
The action is Galilean invariant by virtue of the relation 



A. Galilean invariance 

In two dimensions, the equation of motion 0} and thus 
the action (|18J) have the important physical property of 
invariance under a Galilean transformation to a refer- 
ence frame moving with constant velocity u [Tjj \l$L |20| . 
This property depends on the assumption that the stir- 
ring force is uncorrelated in time, since otherwise there 
would exist a link between a point in space at one time 
and a "corresponding" point in space at a different time. 
Specifically, the action l|18fl is invariant under the trans- 
formation 

ip(t,x) — ► ip(t,x. + ui) +CijXiUj, 
P(t,x) ->p(i,x + ui), 

which induces the familiar transformation of the velocity, 

v(i, x) v(t, x + ui) - u. (44) 

The original Navier-Stokes equation is not Galilean 
invariant because of the friction term, which introduces a 
preferred state of rest; but the differentiation in deriving 
the stream-function equation 0} eliminates the constant 
shift in v H3- 

We would not expect Galilean invariance in d = 2 to be 
preserved by our renormalization method unless the the- 
ory remains Galilean invariant when dimensionally reg- 
ulated. We now show that our formal stream-function 
representation in arbitrary d is invariant under the trans- 
formation H43f) . provided that u lies in the physical sub- 
space. For convenience we use the infinitesimal form 

Sip = tUiViip + CijXiUj, 5p — tuiS/ip. (45) 

In the corresponding variation of the action l|18l) . terms 
proportional to t vanish automatically because they cor- 
respond to a simple spatial translation. The interesting 
terms are those where the t is differentiated (—ipW 2 ip) 
and where ip is varied by eijXiUj. Neither of these affects 



u • k ^jITV, k)=mx VT^>, k; 0, k')|*=o; (48) 

both sides equal u • k k 2 . The same relation must then 
hold between the exact 1PI functions: The coefficient of 
k 2 uj in T^ p equals the coefficient of — i(kf — fcf)ki x k 2 
in Y^^ p . Since the former is not renormalized, neither is 
the latter. 



B. Fluctuation-dissipation theorem 

For zero friction (p = 0), the action (|18|l is equivalent 
to model A of FNS [lj , who note that it obeys detailed 
balance and thus is subject to a classical fluctuation- 
dissipation theorem (FDT) ■ A complicated diagram- 
matic argument demonstrates that the FDT is preserved 
to all orders of renormalization [22^. Here we reach this 
conclusion by obtaining the FDT from an exact symme- 
try of the action in arbitrary d. Under the formal discrete 
transformation 

p^p- 2iip, (49) 

the action with p — changes only by reversing the sign 
of the viscous term. The change in the interaction term 
vanishes upon integration by parts, just as in deriving 
conservation of energy. To restore the sign of the viscous 
term, we further perform a complete time reversal, 

t -> -i, ip -tp, (50) 

which naturally reverses the sign of the dissipation. The 
net effect is the transformation 

p{t) -» p(-t) + 2iV(-t), iKt)--V(-t), (51) 

a generalized time reversal that is its own inverse and 
leaves the action invariant. 

The FDT is derived from this symmetry by expressing 
the invariance of (pip): 

(W> =-<#>- 2 W>). (52) 
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By invariance under time translations and spatial rota- 
tions, negating the times in a two-point correlation func- 
tion is equivalent to interchanging the points. As a first 
application of the FDT, we make use of the theorem that 
the equal-time correlation function (pip}= — (ipp}= is ex ~ 
act at tree level j^- Thus, for p = 0, the exact equal- 
time stream-function correlation is 



{•H>)= = 3*«P^>= + (#}= 



<#>o) 



-2ig~ 



2ir uj 2 + g~ 2 k 4 



The energy spectrum, in the units implied by Eq. i|17[l. 
is then 



(54) 



This is an equipartition spectrum, in agreement with 
FNS model A 

For p = 0, the FDT also implies that, if (jnp) 
and thus (tpp) are made finite by renormalizing g, then 
{ipip) is likewise finite, without the need for field or time 
rescalings. As we have seen, the three-point 1PI function 
is automatically finite. Hence the anomalous dimensions 
vanish exactly for p — 0. But in minimal subtraction, 
the counterterms for rescalings and for dimensionless cou- 
plings are independent of the mass parameter |2j. We 
conclude that in minimal subtraction, even with friction, 



7V = 7 P = It = 0, 



(55) 



and /3(g) depends only on g. Indeed, we have seen these 
statements verified to one loop in Sec. Ill Dl 



C. Renormalization-group flows 

We have shown that the anomalous dimensions van- 
ish, and that the renormalized couplings in exactly two 
dimensions obey 



»~r = p{g) = ^r- 

ap oZtt 



0(g 5 



P : 



p V 



(56) 
(57) 



Thus p is very simply related to p and can be used to 
parametrize it. In the RG flow to low wavenumbers, p 
steadily increases as friction becomes more important. 
Meanwhile, g flows in its own characteristic way regard- 
less of the value of p\ the most we can say is that once 
g becomes small, it decreases further and further, ap- 
proaching zero in the IR limit. The solution of Eq. Q 
in this limit is 



1(P) 



16tt 



\n{k g /p) 



(p < kg), 



(58) 



where k g is the scale at which g becomes large. 

In the IR limit, we expect good accuracy from per- 
turbative results such as the tree-level expression for the 
energy spectrum, 



E(k) 



■nk 



8ir(gp + k 2 )' 



(59) 



We might suppose that the true asymptotic behavior is 
given by replacing g with the renormalized coupling g(k). 
This would follow from RG theory if the loop corrections 
to Eq. I)59|l contained \n(k 2 ). But with gp providing an 
IR cutoff, the diagrams are regular as k — > and instead 
contain Irt(gp). Hence, for the purposes of correlation 
functions, the RG flow effectively halts for wavenumbers 
below the "mass" y/gp. If m is the suitably renormalized 
value of this mass, then as k — > we expect that 



E{k) 



k 3 



8irpg(m) 



(60) 



On the other hand, in the bare tree- level result (|59|l . g 
can be interpreted as a coupling renormalized at a very 
high wavenumber (the forcing scale or UV cutoff). With 
/3(g) > 0, we have g(m) < g, and so the RG result lj6TJ|) 
gives a greater E(k) at low k. Whereas the bare tree- 
level calculation ignores all interactions between scales, 
the effect of renormalization is to place more energy and 
dissipation at low k (and therefore less at high k), con- 
sistent with the inverse cascade. 

From Eq. 1)57(1 . any RG fixed point (g*,p*) must have 
p* = 0. Not only do the anomalous dimensions vanish at 
any fixed point, but for p — we have the exact equipar- 
tition spectrum E(k) cx k , whether g is at a fixed point 
or not. It is clear that, despite the suggestive evidence of 
scale invariance of the inverse cascade, an RG fixed point 
in our framework cannot be the explanation for the ob- 
served A: -5 / 3 spectrum. Nevertheless our theory contains 
all the essential ingredients that have produced the sta- 
tionary inverse cascade experimentally and numerically. 
A natural explanation is that the fc -5 / 3 spectrum arises 
from the nonperturbative behavior of correlation func- 
tions at p > and at large values of g that flow rapidly 
with scale. Although it is far from obvious how such an 
RG flow can produce approximate scale invariance with 
an effective anomalous dimension, we are motivated to 
seek hints about the theory's strong-coupling behavior. 
The first step, which can be useful for more mundane 
purposes as well, is to extend the renormalization to the 
next order of perturbation theory. 



IV. TWO-LOOP RENORMALIZATION 
A. Renormalization prescription and diagrams 

Calculating (3(g) consistently to two loops requires a 
precise specification of the renormalization scheme. Re- 
markably, though, as long as (3(g) depends only on g, the 



(3 function to two loops (but no further) is independent 
of the particular scheme chosen 0. As in Sec. Ill L)l wc 
take p > 0, compute two-point 1PI diagrams to 0(q 4 ) 
expanded about e = 0, and renormalize by minimal sub- 
traction. The two-loop expression for the bare coupling 
in terms of the renormalized coupling is 



-g 3 aie~ 



-a'e- 1 



) + 0(f)}. (61) 



With zero anomalous dimensions, (3(g) is the only RG 
function to be determined and we need only compute a 
single 1PI function to two loops. Below we will choose 
— T pp because its two-loop diagrams are technically sim- 
pler than those of -I*. 

We have seen that gp, which has dimension 2 indepen- 
dent of e, acts as the IR cutoff for wavevector integrals. 
By dimensional analysis, T pp to 0(q 4 ) has the form 



t pp = q 4 



r'+gigp)-^ 2 (^ + &i + 0(e)) 

(^ + b -f+0(e^+0(g 5 ) 
r 1 +g(- + b[ - \h \n{gp) + 0(e) 



9 j 



b' 2 - b 2 ln(gp) 



+ 0(e Q ))+0(g 5 ) . 

(62) 





FIG. 3: Nonvanishing two-loop diagrams for the two-point 
1PI function — T pp . Along with the overall divergence, each 
diagram has at most one divergent subdiagram (bold lines). 



B. Analytic calculations 



To simplify the two-loop calculations we formally set p = 
g~ 4 and restore the ln(gp) terms at the end. 

Multiloop diagrams can be characterized by their over- 
all UV divergence (as all loop frequencies and wavcnum- 
bers go to infinity together) and their subdivergences (as 
a subset of loop frequencies and wavenumbers go to infin- 
ity while the rest remain finite) . Overlapping divergences 
occur when two or more divergent subdiagrams share a 
propagator; this is a definite complication in evaluating a 
diagram, though not insurmountable 2]. Some two- loop 
diagrams for the 1PI function — T^ p contain overlapping 
divergences, so we choose to calculate the diagrams for 
— T pp , which fortunately do not (Fig. [3J|. A three-point 
1PI subdiagram is divergent only if its external lines are 
ipipPi t w0 w iggly and one plain; Galilean invariance does 
not eliminate this subdivergence when we treat each two- 
loop diagram separately, because the finiteness of — rf^ p 
results from a sum over distinct one- loop diagrams. In 
Fig. |3 we have omitted one conceivable diagram contain- 
ing the 1PI subdiagram — T^, which vanishes as noted 
in Sec. lITCl 

Previous two-loop calculations of similar complexity 
have been made for different problems using diagrams 
of the same topology: the Burgers-Kardar-Parisi-Zhang 
equation for interface growth |23j and the Navier-Stokes 
equation in more than two physical dimensions, where 
there are fewer divergences |24j . 



We have programmed MATHEMATICA [25| to automate 
the steps in evaluating each two-loop diagram. The loop 
frequencies are integrated one after the other by adding 
the residues of all poles in the upper half-plane. The re- 
sulting integrand, containing the external wavevector q 
and the loop wavevectors ki.2, is expanded to 0(q 4 ) and 
averaged over directions of q in the physical subspace. 
The numerator of the integrand is now rife with the pro- 
jector both from averaging qtqj and from applying 
Eq. J7|) to cross products. To eliminate Ojj, we average 
the integrand over orientations of the physical subspace 
within (i-dimensional space: We introduce an orthonor- 
mal physical basis, write Qij = ciidj + bibj, average over 
directions of a in the d — 1 dimensions perpendicular to 
b, and then average over directions of b in (i-dimensional 
space. The result is a function of kf, k%, and ki • k2 to 
be integrated over d d k\ d d k2- 

Because we need only the divergent parts of the two- 
loop diagrams as indicated in Eq. I|62|l . we eliminate nu- 
merator terms that produce neither an overall divergence 
nor a subdivergence by power counting. We seek to inte- 
grate first over the loop wavevector (say ki) associated 
with the subdivergence (if any). To simplify the denom- 
inator, Feynman parameters |2| are introduced, and ki 
is translated by a multiple of k. 2 - After the numerator is 
isotropized, the ki integration is done analytically; there 
remains an integral over Feynman parameters and over 
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the magnitude k<z . From Eq. I|62(l , the integrand (exclud- 
ing g 3 q 4 ) has dimension —1 — 2e. The part that behaves 
like fc^ 1_2c as ki — > oo is subtracted and separately inte- 
grated over k2, producing another e^ 1 factor; all the e~ 2 
terms arise here and are found analytically, but the e _1 
terms involve intractable Feynman-parameter integrals. 
Meanwhile, the remaining subtracted integral converges 
at k2 = oo but contains analytically integrable e _1 poles 
from the ki integration. 

Adding the divergent parts of the eight diagrams gives 



t^PP 3 4 



1 



2,2 



20487r 2 e 
-21n(47r) + 27 E 



X 



40967r 2 e 



(63) 



Here X is a sum of Feynman-parameter integrals of com- 
plicated rational functions with integer coefficients; thus 
we expect that X may be a rational number. Though we 
are unable to calculate X analytically, Eq. (|63[) already 
displays some important features: If the e -2 term were 
different, the renormalization performed below would be- 
come inconsistent [23; and if the coefficients of 7e and 
ln(47r) were different, these constants would (contrary to 
expectation) appear in the two-loop /3 function. 



C. Numerical calculations 

We have evaluated X numerically by multidimensional 
Monte Carlo integration, using the technique of impor- 
tance sampling 26] to select more points in the "corners" 
of Feynman-parameter space (with one parameter near 
1 and the others near 0). This technique improves the 
statistics because the integrands tend to diverge in these 
corners (but slowly enough that the integrals converge). 
We treat the integrals that make up X separately, since 
they vary in complexity and in number of Feynman pa- 
rameters. To optimize the precision of the result for X in 
a given computation time T, we reason as follows. The 
contribution to X from diagram i, evaluated with rii sam- 

— 1/2 

pie points, is obtained with a precision Si — Oin i in a 
time ti = TiUi, where Ui and are characteristics of the 
integrand and the computer. Constrained optimization 
shows that we achieve the best overall precision 



in the time 



by choosing 



r = J> 



1/2 ^ 1/2' 



(64) 



(65) 



(66) 



Short runs are made to estimate <7, and r,-, and then the 
high-precision integrations are performed with this plan. 
Our computations for a total of 3.6 x 10 8 sample points 
yield 



X = -3.995 ±0.005, 
strongly suggesting that the exact value is 
X = -4. 



(67) 



(68) 



With the one-loop result l(3T)l and the two-loop result 
we have 



ppp _ qi 



327T6 



ln(ffp/47r) + 7e 

647T 



0(e) 



1 



2048ttV 



21n(gp/47r) + 27E + l , . 
40967r 2 e 1 ' 



0(g 5 ) 



(69) 



By substituting Eq. (|6f I) and requiring a finite expression 
at e = 0, we determine 



at 



1 

32^' 



f 



«2 



2048tt 2 : 



4096tt 2 ' 



(70) 



The RG invariance of the bare coupling g finally gives 



/?(5) = -k9 + Sr 



9" 



32tt 2048tt 2 



+ 0(f). (71) 



This is reminiscent of the (3 function in four-dimensional 
4> theory, where similarly the one-loop term is positive 
and the two- loop term is negative • A positive func- 
tion that grows too quickly (faster than linearly) at large 
coupling raises the question of whether the coupling be- 
comes infinite at a large but finite renormalization scale. 
In our theory this would suggest an absolute limit on the 
extent of the scaling range for any inverse cascade. To 
avoid such a fate, the expansion of the f3 function must 
contain many negative terms to slow its initial supcrlin- 
ear growth. It is pleasing to find such a term already at 
two loops. 



V. INVERSE-CASCADE RANGE 
A. Inverse-cascade phenomenology 

The initial prediction of the two-dimensional inverse 
energy cascade Q assumed zero friction and small but 
nonzero viscosity. Kinetic energy, continually injected at 
high wavenumbers, is expected to cascade down through 
a quasisteady inertial range extending to lower and lower 
wavenumbers as time passes. Within this range, the as- 
sumption of scale invariance implies the energy spectrum 



E(k) = C£ 2 / 3 k~ 5 / 3 , 



(72) 
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where C is the two-dimensional Kolmogorov constant 
and £ is the rate of energy injection per unit mass. Nu- 
merical simulations with zero friction 0, 0, confirm 
this spectrum until the cascade approaches the minimum 
wavenumber associated with a finite system. Both non- 
stationary behavior and finite-size effects, however, lie 
outside our theoretical framework, which treats a fluid of 
infinite size that has been stirred for an infinite time. In 
the absence of friction, such a fluid has an equipartition 
spectrum as shown in Sec. 1111 Bl 

When friction is introduced, it is natural to expect a 
mirror image of the viscous energy dissipation at high 
wavenumbers in the three-dimensional direct cascade: A 
stationary inverse cascade should develop, with dissipa- 
tion by friction at low wavenumbers k ~ kf r and with the 
spectrum <|72[1 at wavenumbers k 3> fcfr where dissipation 
is unimportant. Dimensional analysis gives 

kfr = f-^V/ 2 , (73) 

where a is the friction coefficient 0, . Several numer- 
ical simulations @, @ and laboratory experiments 0, |2j| 
confirm this picture. Other numerical studies obtain sim- 
ilar results but are more difficult to relate to our frame- 
work, since they modify the friction term by removing 
derivatives |7j[ or by applying friction only below a cutoff 
wavenumber |30j ] . We are less concerned about the com- 
mon practice of adding derivatives to the viscous term 
(hyperviscosity), because such a modified term is irrel- 
evant and the RG flow should introduce a normal vis- 
cous term to replace it. But friction modified with in- 
verse derivatives (hypofriction) is certainly relevant and 
is believed to alter the dynamics of the inverse cascade 
0,0. In an extreme case, hypofriction with eight inverse 
Laplacians destroys the /c~ 5 / 3 spectrum |31|. Hypofric- 
tion is intended to confine dissipation explicitly to low 
wavenumbers; ordinary linear friction accomplishes the 
same thing more gently, but makes it difficult in practice 
to achieve an inertial range |2£j. We have used linear 
friction because it is physically realistic and leads to a 
local field theory. 

For convenience in relating the observed inverse cas- 
cade to our field theory with the action (JTSJ, we adopt 
units of time in which 

gv = Dl'\-^ = l, (74) 

so that the rescalings in Eq. (|17|) are trivial. Then the 
kinematic viscosity is 

v = g-\ (75) 
the friction coefficient is 

a = p, (76) 
and the force correlation is 



In contrast to the formal methods of dimensional reg- 
ularization and minimal subtraction, a simple UV cutoff 
renders our theory finite in a physically meaningful way 
while still preserving the symmetries noted in Sec. II I II 
A cutoff would have been inconvenient for the two-loop 
calculations of Sec. I1VI but it is appropriate for under- 
standing experimental and numerical results on the in- 
verse cascade. The local, renormalizable action (I18fl ap- 
plies only with a UV cutoff A below the wavenumbers 
of the external forcing 0] • A corresponding renormal- 
ization prescription is obtained by staying in exactly two 
spatial dimensions and writing the coupling g in Eq. i|18|) 
as a cutoff-dependent quantity g(A), such that the long- 
distance behavior is independent of A. No cutoff depen- 
dence is needed for p, because the friction term is not 
renormalized, as we saw in Sec. Ill Dl 

Special properties of minimal subtraction allowed us 
to conclude (Sec. IIII Bl) that in that scheme the anoma- 
lous dimensions vanish and (3(g) depends only on g. We 
may expect that in the cutoff scheme these statements 
remain approximately true, at least for small values of 
the dimensionless friction parameter A~ 2 p. Because the 
fi function is scheme independent through two loops 0, 
we have 

This result is directly useful for weak coupling, but can 
only hint at the possible strong-coupling behavior. We 
will now show that a specific strong-coupling form of the 
(3 function is required for a fully developed inverse cas- 
cade, and that it can be naturally interpolated with the 
two- loop result (|78ll . 

B. Strong-coupling behavior 

The key condition for the ideal inverse cascade is that 
the dissipation of energy is dominated by friction and is 
almost totally confined to low wavenumbers. This means 
that the dissipation rate 

/■A pA t 

£ = 2p dkE(k) = 2p dk-—^k 2 M} = (79) 
Jo Jo ( 27r ) 

should be independent of the UV cutoff A, because a 
change in A produces neither a renormalization of p nor 
a rescaling of -0. Under stationary conditions, £ equals 
the rate of energy injection, given by |l4j 

f A irk A 4 
Jo ( 27r ) 16tt3(A) 

in terms of the force correlation (|77|l . For £ to be inde- 
pendent of A, we must have the strong-coupling behavior 



D{k 2 )=g^k 2 . 



(77) 



g(A) ~ const x A 4 , 



(81) 
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corresponding to the asymptotically linear /3 function 

0(g) (£->oo). (82) 

We now attempt to connect this form with the two-loop 
function [frSJl. 

Perturbative expansions such as Eq. H78|) are usually 
divergent, but a Borel transformation is expected to pro- 
duce a finite radius of convergence about zero coupling 
. We see that the true expansion parameter is g 2 , and 
the alternatively normalized action (|15fl makes it clear 
that the theory is unstable for g 2 < 0. Thus we write 



9 Jo 9 

For the perturbation series 



A{g 2 ) = Y J A n c 



the Borel transform is 



*M = £ nr 



— z n . 



In our case, 



32tt 4096tt 2 



(83) 



(84) 



(85) 



(86) 



The Borel transform often has poles on the negative 
real axis associated with instantons [2j, but (3(g) is well- 
defined from Eq. (|83|l as long as B(z) is regular on the 
positive real axis. A simple and suitable rational (Pade) 
approximant to B(z) is 



B(z) 



32ir + yz 



(y>0); 



(87) 



fortunately, the choice y = \ agrees with Eq. (|86() . We 
thus take 



B(z) 



4z 



1287T + Z' 



but we do not here attempt to study the possible instan- 
ton solutions corresponding to the pole at z = — 128-7T. 
Eq. then gives precisely the desired asymptotic be- 
havior 



13(g) -4.9 (ff-foo). 



(89) 



We now ask whether the observed k~ 5 ^ 3 energy spec- 
trum H72[) is consistent with our theory, although we are 
unable to derive it systematically. We conjecture that 
strong-coupling effects produce the spectrum 



for 



E(k) ~ £ 2 / 3 fc- 5 / 3 

k>ktr =£- 1/2 p 3/2 . 



(90) 



(91) 



From Eq. H80(l . the running coupling is 
5 (A) -f^A 4 . 



(92) 



For nonperturbative effects to be operative down to the 
wavenumber kf T , we must have 



(93) 



In the borderline case where g(kf r ) ~ 1, we have fcf r — 
5 1 / 4 p 1 / 2 , and we can match Eq. H90(l in order of mag- 
nitude with the perturbative energy spectrum (|59(l : 



E(k b ) - £ 2 / 3 fc~ 5/3 ~ £ 1 ^ 



k 3 



p + K 



(94) 



Even when g(kf r ) ^> 1, we may guess from Eq. lj5*9"|) that 
in order of magnitude 



E(kfr) 



k 3 
A "fr 



P5(fcfr) + fc 2 



£ 3 / 2 p- 5/2 , (95) 



and this again matches Eq. I|9UII at fcf r . 

In the case 5(fcf r ) 3> 1, the running coupling eventually 
becomes ~ 1 at a lower wavenumber 



kg — £ ^ , 



(96) 



below which perturbation theory is applicable and the 
energy spectrum is given roughly by Eq. (|60|l . Thus we 
envision a varied but continuous behavior of the energy 
spectrum in the different regions: 

(p^k 3 (k <£!/*), 
E(k) - I £p- 1 k~ 1 (£^ A <k< £-^ 2 p 3 / 2 ), (97) 

|^2/3 fc -5/3 ( k > £ -l/2 p 3/2y 

As a further check, we note that the contribution to the 
energy dissipation rate l|79[l for each of the three regions 
is — £ (up to a logarithmic factor for the middle region) . 
This suggests that a fully developed stationary inverse 
cascade has several distinct dissipation ranges. Testing 
by experiments and simulations is not straightforward, 
because finite-size effects may become important before 
the lower- wavenumber ranges are reached. 



C. Energy flux and third moment 

We have described the ideal inverse cascade in terms 
of external forcing confined to high wavenumbers (as- 
sumed in Sec. IIIBl! and energy dissipation practically 
confined to low wavenumbers (made plausible in Sees. 
IIIICI and IV B|) . This separation is equivalent to the con- 
ventional criterion of an inertial range with a constant 
energy flux. However, unlike the energy spectrum (a sim- 
ple correlation function) , the energy flux in wavenumber 
space is not invariant under the RG flow and is not given 
straightforwardly by our local field theory. Renormaliza- 
tion introduces the forcing (|77|l . which is nonzero even 
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for wavenumbers in the inertial range and so makes the 
flux appear nonconstant. This effective forcing simply 
represents the energy transfer from wavenumbers k > A 
that have been integrated out. 

We can resolve the flux problem by working instead 
in physical space, where the effective forcing is a differ- 
entiated delta function that is zero at finite distances, 
including the inertial range of lengths. Thus we inquire 
whether the physical-space energy flux 

e(r) = -|V r • (|<Mr)| 2 <Mr)) (98) 

is preserved under renormalization. The flux is given in 
terms of the third moment of the velocity increment 

<Jv(r) = v(x + r) - v(x) (99) 

and is independent of x by homogeneity. Under station- 
ary conditions, the Navier-Stokes equation Q yields [33] 

e(r) = {vV 2 r - a)(v(x) ■ v(x + r)) + ±C(r) 

= [5(A)- 1 V?-/ 3 ]<v(x)-v(x + r)), (100) 

where C(r) is the vanishing physical-space force correla- 
tion, and we have used Eqs. (|?5jl and J7BJ). 

The right-hand side of Eq. (|100fl can be interpreted 
as minus the rate of energy dissipation at length scales 
larger than r. Since (v(x) • v(x + r)) is an RG-invariant 
correlation function, only the viscosity g(A)^ 1 introduces 
cutoff dependence. In accordance with our argument in 
Sec. IV Bl as long as the energy dissipation is dominated 
by friction, the dissipation rate is independent of A, giv- 
ing a well-defined energy flux 

e( r ) = -P< v ( x ) ' v ( x + r )) (101) 

proportional to the Fourier transform of the energy spec- 
trum. And if this spectrum is almost entirely concen- 
trated at low wavenumbers k < kf r , it follows that in the 
inertial range (r <C k^ 1 ) the flux is constant, 

e(r) ~ e(0) = -p(v 2 ) = -£. (102) 

D. Inter mittency 

A striking feature observed in the inverse cascade is 
the approximate scale invariance of inertial-range veloc- 
ity correlations (structure functions). The k~ 5 ^ 3 energy 
spectrum 1(72(1 gives for the second moment Q 

(|5v(r)| 2 > cx£ 2/3 r 2/3 , (103) 

and the constant energy flux (|98|l gives for the third mo- 
ment 

(|5v(r)| 2 <Sv(r)) cx£r. (104) 

In our field theory, since the anomalous dimensions van- 
ish, such moments can be written in the form 

((6v) n } = r- n f n (g(r- 1 ),pr 2 ), (105) 



where r~" is the kinematic scaling and /„ is a function 
of the dimensionless running couplings. Eq. I|105|) should 
reduce to the observed forms in the inertial range 

r<k£ 1 =£ 1 ' a p-*' a , (106) 

subject to the condition for a fully developed cascade, 

g(k h ) ~ (£-V) 3 » 1. (107) 

Conditions (|106f) and (|107f) are equivalent to 

(pr 2 )- 2 « gir- 1 ) « (pr 2 )" 3 . (108) 

To achieve complete scale invariance of inertial-range 
structure functions (absence of intermittency), we must 
have 

/nC^r-^pr^ocflCr- 1 )-/ 3 (109) 

for the range of arguments ((108(1 : then f n oc £™/ 3 r 4 ™/ 3 5 
and 

({Sv) n ) cx£ n/3 r"/ 3 . (110) 

Eq. 1109(1 requires a specific nonperturbative behavior of 
velocity correlations at strong coupling, mimicking the 
effect of an anomalous dimension. RG theory alone does 
not constrain the functions f n . The well-established sec- 
ond moment (|103|l and third moment (| 104(1 strongly sug- 
gest that an exact calculation in our theory would yield 
the required behavior of ji and f^. But in the absence 
of a field-theoretic reason why Eq. 1(1091) should persist 
for n > 4, high-order structure functions may generically 
be expected to violate scale invariance and produce in- 
termittency. In sum, we would not be at all surprised 
by observations of intermittency in the inverse cascade, 
but a seeming total absence of intermittency would raise 
questions about unknown properties of our theory that 
enforce effective scale invariance. 

Some numerical studies of higher moments 0, deal 
with a nonstationary inverse cascade and find no signs 
of intermittency. Such results are not directly relevant 
to this paper but have been explained theoretically [(| 
based on the growth of the inertial range with time. For 
the stationary inverse cascade, laboratory experiments 
0, |3^| reveal no evidence of intermittency. The results 
of stationary numerical simulations, however, are mixed: 
A study using linear friction |8| confirms the absence of 
intermittency, but a study using hypofriction with one 
inverse Laplacian Q obtains intermittency that is de- 
scribed as strong and as similar to that observed in the 
three-dimensional direct cascade. Surprisingly, results 
from the latter simulation are also presented in a subse- 
quent paper |33J where the intermittency is described as 
insignificant. The evidence on intermittency in the sta- 
tionary inverse cascade is unclear, and further numerical 
studies would be useful to resolve the question. Results 
exhibiting strong intermittency are most natural, from 
the viewpoint of our field theory. 
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The mild hypo-friction used in the numerical study ob- 
taining strong intermittency Q is unlikely to alter the 
qualitative structure of our theory. While the inverse 
Laplacian renders the initial Navier-Stokes equation JTJ 
nonlocal, the differentiated form and the field-theory 
action Ijl8|) are still local in terms of the stream func- 
tion ip. The hypofriction term certainly violates Galilean 
invariance, but this symmetry holds for zero hypofric- 
tion and so its implications for the RG functions persist 
in minimal subtraction even when hypofriction is added. 
The same arguments used for linear friction in this pa- 
per suggest that intermittency should be expected with 
hypofriction as well. It would be interesting to perform 
a direct numerical study of the effect of modified friction 
on intermittency in the stationary inverse cascade. 

VI. GENERALIZED BURGERS EQUATION 
A. Dimensional continuation 

As an example of an alternative statistical fluid model 
to which our RG methods can be applied but which ex- 
hibits different behavior, let us briefly consider the one- 
dimensional Burgers equation [34[ 

v + vVv - vV 2 v + av = f. (Ill) 

Here v is the velocity field of a fluid without pressure, / 
is the force per unit mass, v is the kinematic viscosity, 
and a is the friction coefficient (not normally included in 
the Burgers equation but useful in controlling the long- 
distance behavior). As with the two-dimensional incom- 
pressible fluid, we assume that the forcing is confined to 
a band of high wavenumbers, and we study the response 
at lower wavenumbers. 

FNS found that the UV-stirred Burgers equation, 
like the Navier-Stokes equation, has a dimensionless cou- 
pling in two spatial dimensions and can be analyzed by 
an e-expansion in d = 2 — e. But their continuation of the 
Burgers equation to d > 1 does not preserve important 
one-dimensional properties, and their e-expansion is not 
fully consistent [ll|3|. With v = a = f = 0, Eq. JTTH> 
yields conservation of the "energy" 




(112) 



This E is not proportional to the physical energy of a 
pressure-free fluid, because it does not account for the 
varying density; nevertheless, E is conserved in = 1, 
and leads to a fluctuation-dissipation theorem (FDT). It 
is not easy, however, to generalize Eq. (|lllfl to d > 1 so 
that a similar "energy" is conserved, while maintaining 
other key properties such as Galilean invariance. This is 
the task we now address. 

If we neglect dissipation and forcing, the conservation 
of E in d = 1 follows from the relation 

= vv + iV(u 3 ) = v(v + vVv). (113) 



The analogue in d > 1 that would yield conservation of 
E= \J d d x|v| 2 (114) 

is 

= v • v + AV • (w 2 v) 
= v [v + A{V ■ v)v + (A - B)V(v 2 ) + 2£(v • V)v] 
= v w. (115) 

Unfortunately, whatever the choice of the constants A 
and B, the quantity w is not Galilean covariant and is 
not a suitable generalization of (v + vVv). 

To remedy this problem, we impose the potential-flow 
condition 

v = V?/>, (116) 

as is commonplace when considering a multidimensional 
Burgers equation. Then, because 

E = J d d xw w = J d d xi){-V -w), (117) 

a scalar equation of motion can conserve 

E= \ J d d xViipViil). (118) 

We take 

= -V • w = -V 2 ^ - AV 2 ip\7 2 iP - 2AV i V j ipViV j il) 
- 3AV. 1 ^V. i V 2 V'. (119) 

If A = -|, this expression is covariant under the Galilean 
transformation 

i>(t,x) -» t()(t,x + ut) -u-x, (120) 

which induces 

v(t, x) -> v(i, x + ut) - u. (121) 

Upon restoring dissipation and forcing, we therefore 
propose 

-v 2 v> - v^v^V - iv^vV 

- § Vj VjVViVjV' + - aV 2 ^ = V (122) 

as a fully satisfactory multidimensional Burgers equation 
that reduces to Eq. 1|) in d = 1. Eq. Ijl22|l is formally 
very similar to the incompressible stream-function equa- 
tion I0J , and we will follow our previous analysis closely. 
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B. One-loop renormalization 

As in Sec. [HI we assume Gaussian forcing and intro- 
duce a path integral with the action 



S = J dtd d x[\{-\7 2 p)D{-\7 2 )p 



= J dtd d x[\{-\7 2 p) D(-X7 2 )p + ii/pVV - iapW 2 ^ 

(123) 

where we have again integrated the pipip term by parts. 
Analysis of dimensions in d = 2 — e proceeds as before, 
since all terms contain the same numbers of derivatives 
as for the incompressible fluid. The final Burgers action, 
to be compared with Eq. l(T%|) . is 

S = J dtd d x[\g- l \7 2 p\7 2 p + ig^pV 4 ^ - ippS7 2 i/j 

- ipV 2 ^ - iViVjP (±6ij VfcVVfeV- + IViipVjip)] . 

(124) 

The scaling dimensions are again given by Eq. (|20|l . 

The only change to the Feynman rules in Fig. ^ is the 
vertex factor, which now becomes 



~ L o 



iilkr 



-karVka- 



•fi^+ki-kaj^l+kx-ka) 



(125) 

Thanks to conservation of the "energy" E, the FDT is 



obtained just as in Sec. lIIIBl and together with Galilean 
invariance it implies that the anomalous dimensions van- 
ish. Thus, to find the new one-loop (3 function, we need 
only recompute one of the diagrams in Fig. |2 We calcu- 
late 



rr = igq 



1 



327T6 



Hgp/iir) +7e + | . . 
64tt l£ ' 



(126) 

to be compared with Eq. 1(30(1 . Curiously, the 0(e x ) 
term is exactly the same and again leads to the one-loop 
(3 function 



32tt 



o{f). 



(127) 



But the O(e ) term in Eq. ^ 126(1 . which would enter a cal- 
culation of the two-loop (3 function, differs from Eq. I(30|l . 
and we have no reason to expect an identity between the 
(3 functions beyond one loop. 

Our multidimensional Burgers equation 1(1221) appears 
not to have been investigated previously, and it is pos- 
sible that numerical simulations in two or three dimen- 
sions could reveal interesting and unexpected behavior. 
The focus here, however, is on the extrapolation to one 



dimension (e = 1), where we seek to explain the known 
behavior |3J| of the ordinary Burgers equation 1|111|) . If 
e is positive and small, then the f3 function ((127(1 has a 
nontrivial IR-stable fixed point at 



167re + 0(e 2 



(128) 



By contrast, when a simpler continuation of the Burgers 
equation was used, the fixed poin t appeared to be IR- 
unstable near two dimensions 12, 34]. From our results 
it is natural to expect an IR-stable strong-coupling fixed 
point in one dimension (e = 1). 



C. Comparison with inverse-cascade model 

Due to the nontrivial fixed point, the response of the 
one-dimensional Burgers equation to UV forcing is dif- 
ferent from that of the two-dimensional Navier-Stokes 
equation — even in the absence of friction. At wavenum- 
bers low enough that the fixed point is reached, Burgers 
correlation functions (at equal or unequal times) exhibit 
purely kinematic scaling, since all anomalous dimensions 
vanish. For example, from Eq. ((20(1 . the scaling of time 
and frequency is given by 



dk = -2 - 



k 3 ' 2 . 



(129) 



Meanwhile, in both models, the FDT implies that the en- 
ergy spectrum (an equal-time correlation function) obeys 
equipartition. These conclusions agree with previous an- 
alytic 01 an d numerical [34| studies of the frictionless 
one-dimensional Burgers equation. 

We have argued that the inclusion of friction dramati- 
cally alters the behavior of the two-dimensional Navier- 
Stokes model, from equipartition to a stationary inverse 
cascade. This is plausible only in the presence of a large, 
rapidly flowing coupling g. We noted in Sec. IVlhow this 
RG flow could give rise to a constant energy flux and a 
^,-5/3 spectrum. No such profound effect is expected for 
friction in the one-dimensional Burgers equation. Corre- 
lation functions will be modified at very low wavenum- 
bers < /O 2 / 3 , but as long as g remains close to the fixed 
point, higher wavenumbers will retain the equipartition 
spectrum. 

An open question concerns the strong-coupling behav- 
ior of our generalized Burgers equation 1(122(1 in two di- 
mensions with friction. If the (3 function happens to grow 
asymptotically in the same way as the Navier-Stokes one, 
then the arguments of Sec. 0can be repeated and it is 
at least possible that this Burgers model could exhibit 
an inverse energy cascade. If true, this should also be 
evident if the equation is studied numerically. It is un- 
clear, however, whether this multidimensional Galilean- 
invariant scalar field theory has a direct physical appli- 
cation. 
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VII. DISCUSSION 

We have presented a statistical field theory capable of 
describing the stationary inverse energy cascade in two- 
dimensional incompressible turbulence, and computed 
the RG functions through two loops. By contrast with 
previous RG studies of turbulence, we have taken ad- 
vantage of the nature of the inverse cascade, with exter- 
nal forcing confined to high wave numbers, to work with 
a conventional local field theory. The consequences of 
Galilean invariance and the fluctuation-dissipation theo- 
rem have been systematically derived based on the un- 
derlying symmetries of the action. After taking these 
symmetries into account, we have evaluated the neces- 
sary two-loop diagrams in dimensional rcgularization to 
obtain the two-loop j3 function 

—3 —5 

0(5) = - 9 —^ + 0{g 7 ), (130) 

f KB) 32tt 2048tt 2 vy h v ' 

which is independent of the renormalization scheme to 
precisely this order. 

Because the anomalous dimensions vanish identically 
in minimal subtraction, no RG fixed point can yield the 
observed fc -5 / 3 energy spectrum of the inverse cascade. 
Instead, we have found that the inverse cascade could 
plausibly arise from nonperturbative strong-coupling ef- 
fects in the presence of friction. The apparent anomalous 
dimension of the velocity must arise from the rapid RG 
flow of the coupling. Cutoff independence of the energy 
dissipation rate (or constancy of the energy flux) requires 
the strong-coupling behavior 

0(g) -*9 (S-oo), (131) 

and we have also obtained this in a heuristic way from 
a Borcl transformation of the two-loop (3 function. The 
observed energy spectrum in the inertial range has been 



matched with perturbative results at the wavenumbers 
where dissipation becomes important. Incrtial-rangc in- 
termittency (violation of scale invariance) is generically 
expected because the coupling flows rapidly with scale, 
but the evidence on intermittency from numerical simu- 
lations is mixed. On the other hand, a similar RG analy- 
sis of the one-dimensional Burgers equation confirms the 
simpler behavior of that model, controlled by a strong- 
coupling fixed point. 

Our greatest difficulty is that the inverse cascade ap- 
pears to be intrinsically a nonperturbative phenomenon. 
To make quantitative predictions, such as the exponent 
of the energy spectrum or the value of the Kolmogorov 
constant, it may be useful to combine a high-order per- 
turbative calculation with an appropriate resummation 
method, as in our analysis of the strong-coupling (3 func- 
tion. Of course evaluating additional diagrams with two 
or more loops will be very challenging, especially if the 
finite parts are needed. At the present stage, the most 
intriguing prediction of our theory is that substantial in- 
termittency should be expected in the stationary inverse 
cascade. It would be helpful to have more robust and 
consistent results from experiments and simulations to 
determine whether this expectation is realized, and if it 
is not, to identify the field-theoretic explanation. 
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